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We investigate the asymptotics of Hankel determinants of the form 



N-l 

det 

j,k=0 



dxuj N {x) FJ \(M 



as TV — > 00 with q and (1 fixed, where fl is an infinite subinterval of R and u>n{x) is a positive weight 
on Q. Such objects are natural analogues of Toeplitz determinants generated by Fisher-Hartwig 
symbols, and arise in random matrix theory in the investigation of certain expectations involving 
random characteristic polynomials. The reduced density matrices of certain one-dimensional systems 
of trapped impenetrable bosons can also be expressed in terms of Hankel determinants of this form. 

We focus on the specific cases of scaled Hermite and Laguerre weights. We compute the asymp- 
totics by using a duality formula expressing the N x N Hankel determinant as a 2(gi + ■ ■ ■ + q m )-fold 
integral, which is valid when each qi is natural. We thus verify, for such q, a recent conjecture of 
Forrester and Frankel derived using a log-gas argument. 

PACS numbers: 02.30.Mv,02.30.Gp,02.90.+p,02.50.Cw 



I. INTRODUCTION 

Consider the multiple integral 

~ ~ M m 

HM,N,m,q(v) '■= / VN(zi)dzi ■■ ■ Lo N (z M )dz M |A M (2)| 2 Y[ n l/'» - z A 2q ' 
■' n ^° i=l»=l 

where 

q := (q 1 , . . .,q m ), M : = (Mli • ■ -iVm), 

ujn(z) is a nonzero and continuous weight function, possibly depending on a parameter N, and 

M 
C 

3i< 



A M (z) := det^zl X ) = J] (z k - Zj ) 

l<j<k<M 



is the Vandermonde determinant. As a notational convenience we also define 

Hm,n '■= Hm,n,q,-(')i 

= / u) N {z 1 )dz 1 ... I uj N {z M )dz M \£\ M (z)Y 



(1) 
(2) 
(3) 

(4) 
(5) 



It is known that M-fold integrals of the form JIJ can be identified with the determinant of an M x M matrix. 
Expanding the Vandermonde determinants in terms of sums over permutations, and simplifying appropriately, we 
find 



M-l 



H M ,N,m,q(p) = M] - det 



j,k=0 



dza(z)z j (z*) k 



where z* denotes the complex conjugate of z, and 

m 

a(z) :=uj n (z) - z\ 2qi . 



(6) 



(7) 



i=l 
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When C R, equation © becomes 

M-l 

H M ,N,m. q {n) = M\ det [aj+fe] , (8) 

j,k=Q 

where 

a„ := / dza(z) z n . (9) 

One says that the determinant JHJ is generated by the function JJJ . 

Since the entries in the determinant (jHJ are of the form oy+k, we have thus identified the multiple integral JQ, 
when fl C R, with a Hankel determinant. Had we instead taken fl to be T, the unit circle in C, then according to © 
the entries of the determinant would be a^-fc, and we would thus obtain a Toeplitz determinant. (It is conventional 
when discussing such Toeplitz determinants to set z = e 10 and to define a n as the integral of a(z) z n with respect to 
dO rather than dz. This merely introduces the nonzero factor i e lS which is easily absorbed into the definition of a(z), 
and so such technicalities are not relevant to our discussion here.) 

When z£T, functions of the form JJJ) are known as Fisher- Hart-wig symbols^(aith.oxigh. we remark that they are 
not the most general examples of Fisher-Hartwig symbols). By extension, we can describe the Hankel determinant 
(JSJ) as being generated by a Fisher-Hartwig symbol which is defined on the real line. 

The asymptotic analysis of Toeplitz determinants generated by Fisher-Hartwig symbols is a fascinating and well 
studied subject (see e.gi and references therein), and rigorous results which describe the large M asymptotic behavior 
of Toeplitz determinants generated by symbols of the form Q are known£. There are a number of important physical 
applications of such determinants (see e.g. 3-4-5 ). It is often the case, as we discuss presently, that the quantity appearing 
in applications is actually the integral Q rather than the determinant directly, and when ft C R we are naturally 
lead to Hankel determinants generated by the symbol (0. We discuss below a number of physical applications in 
which the asymptotics of such Hankel determinants is of interest. A rigorous treatment of these asymptotics is an 
open problem, however a conjectured form for the large M asymptotics has recently been reported by Forrester and 
Frankel in£. Very recently, the in — 1 case of this conjecture has been verified when f2 = R and ojn{x) is a Hcrmite 
weight, by using a Riemann Hilbert approach^. In the present work we verify the conjecture of Forrester and Frankel 
when ujn(z) is either a Hermite or Laguerre weight for any m G N, when each qi G N. 



A. Random matrix theory 

The multiple integral Q has a natural interpretation in random matrix theory. Let us consider the ensemble of 
random matrices with joint eigenvalue probability density function (pdf) given by 

1 M 

Pm,n{x) := |A M (x)| 2 TTwjvfo), (10) 

Hm,n 

and whose eigenvalues lie in 17. When Q = R, concrete examples of such ensembles include the ubiquitous Gaussian 
unitary ensemble (GUE), corresponding to ujn(x) = exp(—2Nx 2 ), as well as more general unitary ensembles (UE), 
corresponding to ujn(x) = exp(— NV(x)) with V(x) an arbitrary polynomial of even degree with positive leading 
coefficient (see e.gi£). When Q = (0, oo) an important example is the Laguerre unitary ensemble (LUE), corresponding 
to ojn(x) = x a exp(— 4:Nx), which includes Wishart matrices and the Chiral GUE as special cases (the latter after 
a straightforward change of variables; see e.g£). Setting x = e %e and ojn(x) = 1 in I|1U|) we obtain the joint pdf for 
the eigenphases 9i S [0, 2ir] of the ensemble of random unitary matrices with Haar measure, often called the circular 
unitary ensemble (CUE). For the purpose of computing expectations, the CUE is equivalent to (110ft with CI = T and 
ojn(z) — V* z - 

If we denote the characteristic polynomial of the M x M matrix X, with eigenvalues x%, . . . , xm, by 

M 

Z M {Vi) ■= det(/ii I -X) = \\iui - xi), (11) 

i=i 

then Q corresponds to the following expectation involving the absolute value of such characteristic polynomials 

1\{\Zm^)\A =%^, (12) 
\ , / ti M .N 

' Pm,n 
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where the expectation on the left hand side of l|12|) is with respect to the joint eigenvalue pdf (|10|l . The case M = N 
is generally the case of interest. 

From we see that expectations of characteristic polynomials of the form appearing in (|12(1 are characterized by 
a determinant generated by the symbol J7J; when Q C R, it is a Hankel determinant, and when £1 = T it is a Toeplitz 
determinant. 

A sizable literature on the correlations of products and ratios of characteristic polynomials of random matrices 
from various ensembles has emerged in recent years, see e.giSiiiiiliiiSii&liiiiiiASi, and significant progress has been made 
in the understanding of such objects. Such quantities have applications in diverse fields including number theory, 
quantum chaos and many-body quantum mechanics. These works consider either exact algebraic relations that are 
valid for finite N, or the large N asymptotics in the usual universal microscopic scaling limits. We shall be interested 
not in scaling limits, but in the limit of large N with /i fixed. Investigations of objects of the form (|12[1 in this limit 
have been reported in^ii. 



B. Impenetrable bosons 

A compelling physical motivation for investigating multiple integrals of the form arises from a consideration 
of certain one-dimensional many-body systems of impenetrable bosons. By impenetrability we simply mean that we 
require the wavefunction to vanish whenever two bosons occupy the same point in space. Such systems have been 
receiving renewed theoretical interest recently due to the possibility of their experimental realization in the near future 
using ultra-cold systems of atomic bosons confined in elongated traps; see e.gm&^. Systems of impenetrable bosons 
with certain specific boundary conditions are known to have ground state wavefunctions of the form 

! M 

ip(xi, . . .,x M ) = 7t^= U V9n{xi) |Ajvf[/(a;i), . . . %i, ■ ■ ■ ,%m 6 D C R, (13) 

v° M 1=1 

see e.g^SiSi. Specifically, systems with periodic, Dirichlct, or Neumann boundary conditions have wavefunctions of 
this form, as do systems confined in an harmonic well. Of these, the harmonically confined system is perhaps the 
most relevant to current experiments. 

It is worth emphasizing that the introduction of such zero-range infinite-strength interactions establishes a cor- 
respondence between impenetrable bosons and a corresponding system of free fermions, a fact first noted irA and 
this is one of the primary reasons for current experimental interest in such systems. Indeed, for the specific systems 
mentioned above, were it not for the absolute value surrounding the Vandermonde determinant, (|13f) would define 
the wavefunction for a system of free fermions. This fact implies that certain quantities such as the energy spectrum 
and the particle density are identical in the impenetrable boson system and its corresponding free fermion system. 
Quantities which depend on the phase of the wavefunction however will clearly differ significantly between these two 
systems. 

One such quantity, of great significance, is the n-body density matrix, which for a system of N + n particles is 
defined as 

PN+n ( x l>Vl> ■ ■ ■> x n,y n ) ■= y^ 71 ^) Jd^ 1 ' ' ■ f D d £ N ^( Xl > ■■■> X n,£l,-- ■ ,€n)iP*(V1, ■ ■ ■ ,Vn,S,l, ■ ■ -,£n)- (14) 

A key observation is that wavefunctions of the form i|13[l admit the factorization 



I ~^ n N n 

1 p( Xl ,... >Xn ,Z l ,...,Z N ) = J-^ l[^/^i}\A n [f(x 1 ),...J(x n )}\l[l[\f(x i )-f^ l )\ ■ ^(&,... ,6r), (15) 

V N+n 1 = 1 l=\i=\ 

and inserting l|15|l into Ijl4(l yields 

PjV+nO^yi' •■•■> x n, Vn) = ^ n V9N{Xi) g N {yi)\An[f(xi), f(x n )] A„ [/(j/l), . . . , /(j/n)]| 



H N , N ,2n,q{f{xx), f(x n ), /(j/i), . . . , /(t/n)) 



(16) 



where in the definition of HN,N.2n.q and iJjv+n,./v we have VL = f(D) and 



g=(l/2,...,l/2) 



u N (z)=g N (f- 1 (z))^(z). (17) 
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For the four specific systems mentioned below (|13fl . this lun(z) given in i|17|) is well defined, nonzero and continuous. 

Hence, from ((HJl we see that p^\_ n is characterized by a determinant generated by the symbol (JTJ). For systems subject 
to periodic boundary conditions this determinant will be a Toeplitz determinants^^, whereas for systems confined 
by an harmonic well it will be a Hankel determinant. Indeed, this link between density matrices for impenetrable 
bosons and Toeplitz determinants generated by Fisher-Hartwig symbols was originally one of the key motivations 
for investigating the asymptotics of such Toeplitz determinants 3 ^^. In light of the possible future experimental 
realization of finite one-dimensional harmonically trapped systems of impenetrable bosons, an important theoretical 
question is the behavior of the corresponding density matrices when TV is large and Xi, yi, . . . , x n , y n are fixed. This 
then provides a direct physical motivation for investigating the large N asymptotics of N x N Hankel determinants 
generated by symbols of the form Q . 

Perhaps the most important quantity is the one-body density matrix. The asymptotics of the one-body density 
matrix for a system with periodic boundary conditions can be rigorously established from the asymptotics of the 
corresponding Toeplitz determinant. The leading order behavior of the one-body density matrix in the case of 
harmonic confinement was deduced inSi using log-gas arguments, and has been recovered iri2^ using a rather more 
direct, yet still non-rigorous, approach. The asymptotics of the one-body density matrix in the Dirichlet/Neumann 
case was deduced ini£, again by log-gas arguments, and has now been rigorously proved ir£, by making use of recent 
results in2&2£. 

We conclude our discussion of impenetrable bosons by noting the correspondence between the joint eigenvalue pdf 
(fTUfl and the wavefunction lfT3)l. The correspondence between the joint eigenvalue pdf ifTU)) and the wavefunction for 
a system of free fermions is well known 2 ^. The correspondence between impenetrable bosons and random matrices 
was first noted by Sutherland^, between systems of impenetrable bosons with periodic boundary conditions and the 
CUE; see also^i. The correspondence between impenetrable bosons with Dirichlet or Neumann boundary conditions 
and the Jacobi unitary ensembles (JUE) was discussed inii2£, and a similar interpretation for the LUE was noted 
ir£. Again, the most interesting case from an experimental perspective is the correspondence between the GUE 
and systems of impenetrable bosons confined in an harmonic well, and this particular system has been the focus of 
considerable recent theoretical study, see e.gi^2£*2I and references therein. 



C. Asymptotics of Hankel determinants 



The asymptotics of large Toeplitz and Hankel determinants has been of long standing interest to mathematicians. 
For Toeplitz determinants generated by well behaved symbols, very precise asymptotic results are given by the Szego 
limit theorems (see e.gi^). Toeplitz determinants generated by the symbol Q are not amenable to the Szego limit 
theorems however since {7J has zeros. Inspired in part by applications to impenetrable bosons Lenard 4 *^ (see also 3 -) 
conjectured the asymptotics of Toeplitz determinants generated by symbols of the form Q, and this conjecture was 
subsequently proved by Widom 2 . The asymptotic behavior of Toeplitz determinants generated by J7J), as well as more 
general Fisher-Hartwig symbols, is now well understood (see e.gi). Analogously, the asymptotic behavior of large 
Hankel determinants generated by functions defined on f2 C R has also been the subject of study. This problem was 
addressed by Szego 3 ^ and also Hirschman 34 with a finite interval (see also^). In the context of the UE and LUE 
of random matrix theory, as well as in the context of trapped systems of impenetrable bosons, we are interested in 
case where f2 is infinite, and recently Basor et aS^ have considered the asymptotics of Hankel determinants generated 
by symbols defined on Q — (0,oo). However, a key restriction in these works is that the symbol be nowhere zero, 
and hence they do not apply to determinants generated by J7J. Forrester and Franked have recently conjectured the 
asymptotic behavior of Hankel determinants generated by symbols of the form Q defined on £1 C M. Complete and 
rigorous proofs of their conjectures remains an open problem. As mentioned above, a rigorous proof for the m = 1 
case when Q = R and ujn(x) is a Hermite weight has very recently been reported in&. 

Specifically, Forrester and Frankel 5 consider the behavior of the ratio 

riN, m ,qW ■= —r ) ( 18 ) 



where 



qi + --- + q m , (19) 



as N — ► oo with /i and q fixed. Amongst other results, they consider the case uj?j (x) = e ~ NV ^ with either = R 
or £1 — (0, oo). For a particular choice of such ujn{x) and Q there corresponds the quantity p(x), which, with Pn,n 
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defined as in 1|10|) . equals the limit of 

Pn,n(x, x 2 , ■ ■ ■ , x N ) dx 2 ■ ■ ■ dx N (20) 



as N — > oo with x held fixed; i.e. p(x) is the limiting expected eigenvalue density of the ensemble of random matrices 
defined by Pn,n- We note that p(x) is non-negative and has compact support. In what follows int(supp p) denotes the 
interior of the support of p{x). For a detailed discussion of p(x), and a number of other very interesting alternative 
characterizations of p(x), the reader is referred toi The conjecture reported in 5 (in our notation) is the following. 

Conjecture (Forrester- Frankel). Let m,N € N, q € (— l/2,oo) m and p,i,...,p, m 6 int(supp p), where p(x) is 
defined as above. Furthermore, suppose either Q = R or Q = (0, oo), and take ujn(x) — e ~ NV ( x ) where V(x) is a 
polynomial which is independent of N and which has positive leading coefficient and no zeros in f2. Then 

m m „ 2 / _i_ -i \ 

i—1 i—1 * ' 



(21) 



l<j<fc<m i—1 



Here and in the sequel G refers to Barnes' G-function^ 6 .. 

Despite first appearances, the structure of Ij21|l is actually quite simple. Note that the only ensemble dependent 
quantities on the right hand side are p(x) and ujn(x), and that the dependence oiTi-M^^if) on these two quantities is 
universal. The other quantities on the right hand side are truly universal. We note that the factor G 2 (qi + l)/G(2qi + l) 
occurs also in the asymptotics of Toeplitz determinants generated by Fisher-Hartwig symbols^, and has been discussed 
in the context of moments of random characteristic polynomials^*^. We should note that the actual conjecture 
reported in£ is slightly more general than l|21|) . but (|21|l is sufficient for our purposes. 

As an aside, we remark that the ratio appearing in (|16|l is precisely of the form l|18|l and so the asymptotic behavior 
of all the n-body density matrices for a system of harmonically confined impenetrable bosons follows directly from 

CD- 

The present article focuses on two classical cases already mentioned, the case of the Hermite weight, corresponding 
to the GUE, and the case of the Laguerre weight, corresponding to the LUE. In the Hermite case 

w N {x) = e- 2Nx \ n = R, p( M ) = -yi _ ^, supp(p) = [-l,l] J (22) 

and in the Laguerre case 

co N (x)=x a e~ 4Nx , tt=(0,oc), p ^) = l^l-i t BuppOO = [0,l]. (23) 

The conjecture reported ir£ was deduced by considering the specific examples of the Hermite and Laguerre cases, 
to find the general form in terms of p{x) and wjv(x). The asymptotics for these two cases was deduced by using 
the log-gas analogy to conjecture a factorization of (|18fl . computing the asymptotics of each factor when the qi were 
natural, and then conjecturing an analytic continuation to real q^. 

In this work we show how (|21|l can be proved rigorously for the Hermite and Laguerre cases, when each qi is natural, 
by using a duality formula derived from a general result irA By a duality formula we mean an equation identifying 
the A^-fold integral TiN,m,q{n) with a 2\q\ - fold integral. Our results clarify the origin of the factors appearing in 
(|21|l . A similar approach has been used inS to investigate a particular special case when the weight was of Hermite 
type, in the context of trapped impenetrable bosons. 

Section ^ contains a discussion of the duality formula. Due to the similar nature of the Laguerre and Hermite 
cases, they can both be derived simultaneously. Section ITTT1 then discusses the asymptotic analysis of the 2\q\ - fold 
integral obtained from the duality formula by means of the saddle point method. We show how to deduce the general 
form of all terms and explicitly simplify the leading order term, and thus verify the conjecture (|21|l for the Hermite 
and Laguerre cases when each qi is natural. 
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II. DUALITY FORMULA 

Let us define 



h AN) 



F k ,n(X u . ..,Xk) := — det ^.^(Afc), (24) 



A K (Xi,...,X K ) j,fc=i 

orthogonal polynon 

by the following two conditions 



where {^k }kLo are the monic orthogonal polynomials corresponding to u>n(x) and Q; i.e. they are uniquely defined 



u> N (x) dx 7r (A ° {x)tt[ N) (x) = 0, j^k (25) 

a 

nf\x) =x j +0(x j - 1 ). (26) 

For the Hermite and Laguerre cases, the 7r^^ ._ 1 (a;) can be expressed in terms of the standard Hermite and Laguerre 
polynomials found in Szego's classic book££ as follows 



(27) 



w J 2-^ N +i- i y 2 N-( N +i- i y 2 H N+] _ 1 {V2Nx) 1 Hermite, 

Tjv+j-iW - |(_i)JV+j-i(jv + j - l)\(4N)- N -J +1 L^l^Wx), Laguerre. 

According to Brezin and Hikam^, we have the following very useful identity 

If f N K 

F k ,n(Xi,...,X k ) = — / u N {x 1 )dx 1 ... / uj N {x N )dx N A 2 N (x) TT TT(Aj - x t ). (28) 

Mn,n Jn Jn 

If we restrict ourselves to q £ N m and set 



1=1 i=l 



K = 2q x +--- + 2q m = 2\q\, (29) 



we can consider the confluent limit 



lim := lim . . . lim . . . lim . . . lim lim . . . lim , (30) 

* ^T(m)+2q m -^fJ-m V(,„) + 1 — >f>m Mq± + 2q 2 ~ > L*-2 M qi +1^^2 M qi ^t±l >-l — > ^1 



where 



i-1 



r(i):=^2<ft. (31) 



l=i 

Taking the limit (|30|l of both sides of (|28|l . and using the elementary fact that (/li — xi) 2qi = \fJ,i — xi\ 2qi when qi is an 
integer and ^ and x\ are real, we thus obtain 

hm F k ,n{Xi, ■ ■ ■ , X K ) = tz — , (32) 

* tlN,N 

and therefore 

Hw.mM = J lNN lim F K , N (Xi, ...,X K ), qeN m . (33) 

l?N+\q\,N * 

This is the key relation we need to derive the duality formula for TLN,m,q{n)i all that remains is to take the confluent 
limit of Fk,n- 

For later convenience, we set 

^N+j-x( x ) = QN{x)r N+j -i(x). (34) 
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If we insert <|34|) into (|24|l and take the limit (|30|l , then by factoring the Vandermonde determinant we obtain 



2q k 2 qj K 

\\m.F K , N {\ u ...,\ K ) = lim JJ JJ JJ (A T ( fc)+ i fc - A,.^)^) -1 JJ Civ(Aj) 

l<j<fc<m l k = ll j = l j=l 



lim J| A^(A tW+ i, . . . , A r ( i+ i)) det r N+j ^i(\ k ), 



(35) 



l<j<fc<m i— 1 



A' 

hm]J A^(A T(i)+ i, . . . , A T(i+ i)) det r N+j -i(X k ) 

i=l J ' 



(36) 



To compute the remaining limit in l|36[) we can use the following. 

Lemma 1. Let c(A) denote a column vector with r(m 4- 1) entries, then for i = 1, 2, . . . , m 

1 



lim 



lim 



K(i) + 2 qi -*M K(i)+i~ »W ^2g ( (A T (i)+i, . . . , A r (i) + 2gi) 

x G(2% + l)det [c(Ai) . . .c(A r(l)+1 ) . . . c(A r(i)+2 gJ . . . c(A T(m+1) )] 



det 



c(Ai) . . .c(A r(i )) c(pi) —c(fli) . . . j 2q i -l C (i i i) c (AT(i+l) + l) ■ ■ .c(A r ( TO+ i)) 



d^i v d[i, 
where G is Barnes' G-function. 

Proof. This is easily proven by induction using L'Hopital's rule, and recalling the identity 



(37) 



[Jr(l)=G(n + 2). 



i=i 



Applying LemmaQto (|36|l independently for each set {A T (j) +1 , . . . , A T (j) +2gi } with i = 1,2, ... ,m, results in 

d u - 1 



l<j<fc<m J=l l<i<m 

l<j<2| 9 | 



(38) 
□ 

(39) 



In l|39|) the columns of the determinant are ordered such that one starts with i = 1, writes out the 2qi columns 
depending on Zj, and then moves to i = 2 etc. We remark that we have now already obtained one of the two Barnes 
G-function factors that appear in 1)21(1. 

The special property possessed by the Hermite and Laguerre polynomials that allows us to derive a duality formula 
for TiN,m,q{lJ) for the specific weights 1221) and (|23(l is that they can be expressed in terms of contour integrals. Indeed, 
by suitably massaging the standard results in Szego's book2£ we find 



Cj (N)e 2N ^ / dze- 2Nz ^ + Nz2 / 2 z N+ J- 1 Hermite, 



Cj (N) / dze~ 2Nz ^ 



(z + 2) N+a (I 1 



yN+l 



z 2 



2N 1 



TT i2 N +l 

^N+j-iy. i 



NN+j-i 2 2N +i+ a i 



Laguerre, 

Hermite, 
Laguerre, 



(40) 



(41) 



where in the Hermite case the contour C lies along the imaginary axis and is oriented from — ioo to -Moo, and in 
the Laguerre case C is a closed positively oriented contour which encircles the origin but does not contain the point 
z = -2. 
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It is now straightforward to compute the derivatives required in i|39|) from the contour integrals in (|40|l . Defining 

( e 2N & Hcrmite, 



1 Laguerre, 



and recalling the definition (|34|l we obtain 



jZi — 1 



— r N+j - 1 ( f x i ) = c j (N)d h (N) I dz e- NS ^u{z)z u - l [z s + 



(42) 



(43) 



where 



2/ii z — log(z) — — Hermitc, 
2//j z + log(z) — log(z + 2) Laguerre, 



1 



u(z) := < (z + 2Y 



Hermite, 
Laguerre, 



d k (N) := {-2N) h -\ 



(44) 

(45) 
(46) 

and where 5 — ±1 and c = 0, 1/2 in the Hcrmite and Laguerre cases respectively. 

Our task now is to simplify the determinant appearing in (|39|l by using the contour integral (|43J) . This is achieved 
by the following lemma. 

Lemma 2. If q £ N m and 5 = ±1, and we define 
then 



1 5 = +l, 



(47) 



det 

\<li<2qi 
l<z<m 

i<i<2|?| 



dz T(l)+h e ~ A,s ^<^ i .^^(z rW+;i )z^ ( T ) i +ii ^ w+;i +c 



-i r~<5 



•y- 1 



m r(i+l) 



n n2 +1) n n / ^e-^'^^j^^jn^^m---^)) 

i=1 l i=H i=T (i)+l Jc i=l 



r(fe+l) rO'+l) 

- n n n («> 

l<j<k<m l k = T (k)+l ij=T(j')+l 



Proof. We start with the identity 



dct 

l<h<2qi 
l<2<m 
1<J<2|«| 



dz T ^ i ) + i i g h (z T ^) +u , ^j) [/(z r (j) +i . )] 



7-i 



rn 2q 



Yl n / ffli(2r(i)+i i ,j"i)A 2 |g| (/(zi),...,/(z 2 |«|)) , (49) 



i=l Z i = l 



which is valid for arbitrary integrable functions f(z) and gi^z,^). If we apply 149|) to the left hand side (LHS) of 
and use the elementary fact that 



A n (zi + C, Z 2 + C, . . . ,Z„ + c) = A„(zi, Z2, • ■ ■ , Zn), 



we obtain 



m 2q 



m r(i+l) 



lhs of m = n n / dz ^)+u z%\ h n n e-^*^ . a 2M (^, . . .,4,,). 



(50) 



(51) 



i=l 2 i = l 



i=l i i= T(i)+l 



To proceed further we first note the following two useful identities. 



Lemma 3. If f(zi, . . . , ^aigl) is a totally antisymmetric function of each set of variables {z T (j) +1 , . . . , z T ( i+1 )} ; for 
2 = 1,2,..., m 7 then 



n n / dz r(i)+u z %) +u • • • > z 2\ q \) 

m m / \ 

= II r(2a +1) 11 II dz h\ A 2 qt (z T (i )+1 ,...,z T{l+1) ).f(z 1 ,...,z 2 \ q \) (52) 

i=l ^ ' »=1 \i ! =r(i) + l J / 

Proof. By expanding the Vandermonde determinant and then rearranging the order of integrations we see that 

m / r(i+l) \ 

IT I IT / ^« ^2 ?i (2 T (i) + i, Z r ( i+ i) )./(..., Z r (i) + i, Z T (i) + 2 ?i , ■■■ ) 

i=l \li=T(i) + l J J 

m 2qi 

= I1 E II / ^H^^) ^ 1 )' T! ■/(•■■'^(i) + l'---' Z T(i)+2 9l ,■■■): (53) 

1 = 1 <Ji£S 2qi li = l 

and the antisymmetry of / then implies that the right hand side of (|53(l equals 

n e n / ^(Th^,) • /(•■■ I z r(J)+(Ti( i ) ,...,z T( . i)+0 . i( 2 gs ),...) 

= 11 E II / dZ T{i)+U Z t(i)+U- f(---> Z T(i)+l>---> Z T(i)+2qi ,■■■), (54) 
i=lrJieS 2qi U = l 

where the last equality follows by simply relabelling integration variables. The stated result is now immediate. □ 
Lemma 4. With q, 8 and bs, qi (z) as defined in Lemma\^we have 

A 2qi (zf,...,Z2 qi )A 2gi (zi,...,z 2gi ) = Y[bs, qi (z l ) ■ Al q .(z 1 ,...,z 2qt ). (55) 

i=i 

Proof. When 8 = 1 there is nothing to prove, so take 8 = — 1. Then 

A aw (^,...,^jA 3ft (*,..., 2^) = II (Zj "^ fc ~ Zj) . (56) 

l<j<k<2 qi J ' * 

= (-l) 2 <^ [] — ■ AlJ Zll ...,z 2qi ), (57) 



l<j<k<2 qi ^ k 



2q, 



= n^^ • Al qi ( Zl ,...,z 2qi ). (58) 

□ 



1=1 
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Armed with Lemmas [3] and 0] the proof of Lemma [2] follows at once. Applying Lemma|2]to (|51|l results in 

m m I r(i+l) \ 

LHS of m = II T(2 a + T) II II / dz ^ A ^ (MO+i' • • • ' 

i=l ^ ^ ' i-1 \k=T(i) + l JC J 



m r(i+l) 

x Ii II e~ NS{z '^ } u(z h ).A 2M (zt...,z s 2lql ) (59) 



fit ^| (it V 1 / /» 

II, ,, - ii n «(=„» 



11 T(2 qi + 1) J- 1 11 

i=l y H ' i=l J(=r(i)+1 

m 

X J| A 2g4 (4(i) + 1, 4(t+l)) A 2q z (z T (i)+l, • • ■ , «r(i+l)) 
i=l 

r(fe+l) r(j+l) 

x n n n ^ 

l<j<fc<m j fc=T (fc)+l ij=T(i)+l 

Applying Lemma 2] to the right hand side of l|60[l produces the stated result. □ 

Now we substitute (|43f) into 139JI . factor out the constants Cj(N) and di^N) from the determinant, and apply 
Lemma |21 to finally obtain 

Proposition 1. 

TO ^ TO 

n N , m , q (ii) = n ,7n,q J] rf 2g + i)rf2 - + i) n^^ II (m* - m,-) -4 ** • /jv,™,«(a0, (6i) 

i=l k 5» T ^ ^ 9t T ; j=1 l<j<fe<TO 

where 

to r(i+l) 

A: 



to r(i+l) 

/jv, m , g (M) := IT IT dz h e- NS ^^A 2 2qi (z Ti i )+1 ,...,z T{i+1) ) 



1=1 Z;=T(j) + l ' 

to r(i+l) r(fe+l) r(j+l) 

xn n n n n « 4) 

i=lZ i=r (i) + l l<j<k<m l k =r(k) + l l j= T(j) + l 



(62) 



i/ie function g qi [z] is 



and 



<*N : m,q 



g qi (z) : = u(z)bs, qi (z) (63) 
1 Hermite, 

.(z + 2) a (64) 
; t, Laguerre, 



t(to+1) m 2qi 

n ^nn^-^ ^ 



j=l i=lj 4 =l 



JV+|g|,JV 



Proposition ^ is an exact duality formula when q G N m , expressing the AT-fold integral Wjv,m,g(/-t) in terms of the 
2|g|-fold integral lN,m, q {p)- This allows us to compute the large N asymptotics of TLN,m,q(fj) by computing the large 
TV asymptotics of lN,m,q(p), and the latter can be obtained by using the saddle point method. This is the subject of 
Section nm 

The prefactor hj^ t . m ^ q defined in 165fl can be expressed in terms of the Barnes' G-function by using known results for 
the Selberg Integral, see e.g£2£, and the asymptotics can then be obtained from the known asymptotics of Barnes' 
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G- function 38 . We obtain: 



2 -|9| 2 -3| 9 |/2+i:™ 1 2^ -3|g|/2 /yEgg 2g 2 +|g| 2 /2+|g| N G ( N + 2 ) Hermite 

G(7V+ |g| +2) 

Jltf.m,, = <| (-1)' 9 ' 2- 2 I'l+ES! 29. 2 | 9 | ^(a-1,1) M+E2l 2 9 , 2 

G(N + 2) G(N + a + l) G(A^ + 2|g| + 1) 



G(N + 1) G(N + a+\q\ + l) G(N + \q\ + 2) 



i=i 



1 + 



N 



^rE2i(2g 2 -9i) TT(_i)g» 2 2<? ! 2 -2g l7r -2<?, 

i=l 

It will be useful in Sect ion FTTH for us to introduce the notation 



1 + 



1 



N 



where 



h If , m , q = N'Z2*i<to*-«)h 

m 

\l\N 2 -\ q \* JJ 2 2^- 2qi7T -2 qi Hermite5 



h := 



i=i 



[J (-!)'« 2 2 «* 



-2 qi ^-2 qi 



Laguerre. 



Laguerre, 



Hermite, 



Laguerre. 



(66) 



(67) 



(68) 



(69) 



III. ASYMPTOTICS 



Now we begin the task of computing the large N asymptotics of the integral lN,m, q (p) for fixed \x and q. Since the 
only appearance that N makes in (|62ll is in the exponent of e~ N s ( z > tii \ this problem is a natural candidate for the 
saddle point method. 

In both the Hermite and Laguerre cases, the function S(z, jij) has two saddle points, z +> i and its complex conjugate 
z*, i . Explicitly 



It is worth noting that in both cases 



Hi + i \J 1 — fi 2 Hermite, 

; < ■ ~ = \ 

-l + Ul 1 Laguerre. 

Pi 



Im{z+,i} = -p(jm), 



(70) 



(71) 



where p(pi) is as defined in ij^t an d (123; f° r the Hermite and Laguerre cases respectively. Both saddle points are of 
equal importance, since 

Re{S(z + . l , ^)} = Re{S(z* +<i , m)}, (72) 

and we deform the contour C through both of them. 

Let us denote the subset of the contour in neighborhoods of z+j and z*, i by Cl+j and respectively, and the 
complement of the union of these two neighborhoods in C by C s , so that 

c = c s un +ii un_j. (73) 

By deforming C appropriately, the dominant contribution of each integral comes from f2+,j and 0-^, and the standard 
arguments of the saddle point method lead to 



m r(i+l) / 

i N ,mM = n n / 

i=l l t=T (i)+l \ Jn 



dz k + / dz u \ e -NS(z h ,^) . A|,.(z T(i)+1 ,...,Z T(i+1) ) 



n +i , 

m r(i+l) r(fc+l) r(j+l) m 

x ii n aM n n n «-4)+ii^ NRc{Si} -o(e-^), w 

i=l l i=T (i) + l l<j<k<m l k =r(k) + l lf=r(j)+l i=l 
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for suitably small e > 0, where we have defined 

Si := S(z +ti ,fii) = S{z* + i ^i)*. (75) 

We would now like to expand out the 2 |g|-fold composition of the sum of the two integrals appearing in iJHJ. To 
achieve this, we first note that the integrand in (|74|l is totally symmetric in each set of variables {z T ^ + i, . . . , z T (i+i)}i 
for i = 1,2, ... ,m. With this in mind we can then apply the following elementary result. 

Lemma 5. If f(z\, . . . , z T ^ + i, . . . , z T (;+i), . . . , z T (m+i)) is a totally symmetric function of the variables 
{z T ( i)+1 , . . .,Z T (i+i)}, then 



r(»+l) / . f \ 

Yl / dz h + I dz h \ /(zi,...,z T ( i)+1 ,...,z T ( i+1 ),...,z T ( m+1) ) 
i 4=T (i)+i \ J n+,i Ja- tl J 



2 >3i /„ \ r(i)+m , r(i+l) 

2qi 



^2{ n i ) n / dz <«- n / dz <«- /( z i'---' z r( l )+i,---,z T(l+1) ,...,z T(m+1) ). 



(76) 



Proof. All terms in the expansion with integrals over f2+ j can be seen to be equal by swapping the order of the 
integrations, permuting the arguments of / in {^ T (i)+i, • • • , ^ r (i+i)}) an( i then relabelling the integration variables 
appropriately. □ 

Applying Lemma [S] to (|74|l results in 



m 2g< , . T{i)+m r(i+l) 

w^nEn n / n / ^ 

z=l n 4 =0 V l/ i<=T(»)+l " /n +-< Z i=T (i)+n, + l '* 

m / r(i+l) \ r(fc+l) r(j+l) 

n n Ai w ( ZrW+1 ,..., Mi+1) ). n n n 

i=l Uj=T(t)+l / Kj<Km! t =r(l;) + l! J =r(j)+l 



T-j- e -2 ?i JVRe{S*} . 0(e- siV ). (77) 



i=i 



Now let us parameterize the integration variables in (|77|l so that the paths Sl± t i become line segments, (of length 
2r] say), centered at z+j and z+ i respectively and lying along the direction of steepest descent 



_ i z +ti + e l9i ti t , T~(i) + 1 < k < r(i) + rii, 
"' -~ ' z* +>i + e- ie H h , T{i)+n i + l<l i <T{i) + 2q i 



The angles Bi are chosen in the usual way so that with 

._ S"(z +yi , ixj) ^ 2i9i _ S{zj^jm) 2 ,„ 
ai . e e ^(jj 

we have Oj G (0, oo), and so 



, _ s" (z+,i, m) 2 _s"{z* +ti ,m) 



with ti i £ [—1], rj\. With this convention the integrals through z+.j are oriented in the negative direction, and we will 
compensate for this by introducing the explicit factor (— l) ni . Explicitly, since 

S"(z + ,i,^) _ j^p(^)e^ Arcsi ^ + ^ Hermitc, 
2 1 7r^?p(Mi)e _i7r/2 Laguerre, 



we obtain 



7r — Arcsin(//i) 



7T 

4 



Hcrmite, 
Laguerre, 



I 77 P (/•*») Hermite, 
a t = < 2 

[it [i% p(pi) Laguerre. 

Note that > when ^ G int(suppp). 

In making the change of variables (|78|l in (|77|l . it is convenient to introduce the the following definitions: 



% 6i + 2NIm{Si}, 



^ : = E s^-JIT) (fc + 2)! <*' ~ ^ 



m r(i)+rii t(i+1) 
i=lii=T(i) + l 2;=T(i)+ni + l 

GnW:=n n 

i=lii=T(i)+l 
m r(i+l) T(i)+rii 

i>n(*):=n n n (**.-*») 

i=l fc i =r(i)+n i + l Ji=T(t)+l 
r(fc+l) rO'+l) 

^(*):= n n n 

l<j<k<m l k =r(k) + l lj=T(j)+\ 

F n (t) := G n (t)D n (t)H n (t), 

where * on the right hand sides of l|87[). (|88|l and l|89|l refers to the change of variables J7SJ, and n := (m, 
We also define the (unnormalized) integral operator 



Zi=T(z)+r ' 

Armed with these definitions, we sec that since 



r(i)+p n 

E S := II / dtli eT Na%tl * Al_ r+1 (t T{i)+r , . . . ,t T (i) +p ). 



r(i+l) 



^■nt ( z t(%)+i, z T (i)+ ni ) Al q ._ n . (z r (j) + „ i+ i, . . . , Zt(»)+2 9 J H dz^ 

i i=T (i)+l 



r(i+l) 



Ji=r(t)+1 



and 



m t(i+1) m m m r{i+l) 

II II e~ NS( - Zl * 4li) = Y[e~ 2q * NRc{Si} Y[ e 2t( - qi ~ ni)NIm{Si} Y[ T\ e~ Nait *i 

i=ll i= T(i) + l z=l i=l <=lJi=r(t)+l 

we arrive at the following more compact expression for liv,m,g(/-0 



$»(*), 



I S .m,q{p) — JJe 



-2 qi ATRc{S,} 



m 2<?i /O \ 

n e : 

1=1 n i= v 7 
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The factor (— l) ni results from the fact that we traversed the line through z+ i in the negative direction, whereas the 
line through z* + i is traversed in the positive direction. 

To obtain an asymptotic expansion of Zzv,m,g(/"0 from Q94JI we proceed in direct analogy with the one-dimensional 
saddle point method (see e.g^) and introduce the following generalization of $> n (t), 

m T(i)+m r(i+l) 

*n(*,«):=II E[ J] e u h^h), (95) 

»=lJi=r(t)+l i i =r(i)+n i + l 

where u £ K. 2 ' g ' . For convenience we also introduce the function 



so that with 



we have 



Qn(t,u) := $ n (t,u)F n (t), 



ui t = —Na,it l . for each 1 < h < 2\q\ 



Q n (t,{-Naitl}) = $ n {t)F n (t). 



(96) 
(97) 
(98) 



Let us suppose for the present that u is some arbitrary fixed parameter independent of t, and consider the fcth 
degree Taylor polynomial of Q n (t, u) as a function of t € K 2 ' 9 ' 



i(t,u) 



y - 

0<|a|<fe 



dt 



2|9| 



t=0 



t a + 0{t a ) 



|a|=fe+l 



where a £ Z>q and we use the standard notations a! = a\ \ . . . ctaigi 

! and t a = t* 1 . . . t^ 1 and | a\ = a% 



If in l|99|) we now chose u according to l|97|l we obtain 



Q n (t,{-N ai tl})= A 

0<\a\<k 

+ 0{t a ) 



dt? 1 '"dt. 



$ n (t,u)F n (t) 



2|g| 



t° 



{u h =-N ai tf.} 



(99) 



a 2 \ ql . 



(100) 



q| = A-+1 



where we emphasize that the partial derivatives with respect to t on the right hand side of l|100|) are performed with 
u fixed before making the substitution (|97|l . The effect of constructing the Taylor series in this way is that when l|100|) 
is substituted into (|94|l and the integrations are performed, each term corresponding to a given value of \a\ in l|100|) 
will have the same N dependence. To see this we need to consider the asymptotics of the integral operator (|91|l acting 
on a general monomial, which can be deduced simply by scaling N out of the integral. We thus deduce that 



) 2 -M/2 



(101) 



As a direct consequence of (|101|t . we see that with = —aiNtf, for each 1 < h < 2|g|, we have 

m 

H^lXnl+l^ ru a ' =0(AT-E™ 1 9?-E£ I K-« i ) 2 -H/2^ ( 102 ) 
i=l 

for any a[, . . . , a^y G N; i.e. the left hand sides of (|101|) and (|1U2I) have precisely the same asymptotic dependence on 
N. From (|95l) its clear that when computing the partial derivatives of $„(£, u) F n (t) for a given \a\ various powers of 
u will appear, which after the substitution (|97|l will mean we need to calculate quantities of the form 1102fl . However 
(|101fl and l|102f> tell us that all such quantities have the same iV dependence. This would not be the case if we had just 
naively constructed the Taylor polynomial of & n (t)F n (t). We remark that since the symmetry of the integral operator 
(|91|l implies t a is annihilated whenever \a\ is odd, only integer powers of 1/N actually appear in the asymptotic 
expansion of lN,m,q(.n), despite the appearance of \ot\/2 in the exponent in (|102fl . 
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As a result of the expression (|102fl we can see that the dominant term in the asymptotic expansion of l|94|) occurs 
when both rii = qi for all i = 1, . . . , m, and |a| = in the Taylor expansion (|l(J0fl . In general, the coefficient of the 
term which is of order 1 /N k relative to the leading term is composed of all terms for which 



»=i 

If we are interested only in retaining the leading term the preceding arguments imply that 

1 + 



m e -2 qi NRe{Si} 

lN,m,q{v) = Y\_ ' ^0 

i=l 



Nit 



(103) 



(104) 



where the coefficient Iq depends on /i and q but is independent of N. To obtain the explicit form of Ig we first note 
that, since 3> n (0) = 1, we have 



E« lir+p $„(0)F„(0) 



1 



NP 2 / 2 



F n (0)Z p (a t ) + 0(e-^ a > N ) 



where 



r p 

Z p (a % ) := / dPxA 2 p (x) Y[e~ a * x 
Jv I=1 



7T \p/2G(p+2) 



(2P- 1 ) 



P 2 /2 

a. 



(105) 

(106) 
(107) 



The quantity Z v (a,i) is the normalization of the joint eigenvalue pdf of the GUE and can be expressed in terms 
of the Selberg integral; se o 8 i 28 . If we substitute (|100|) into (|94|l and take the leading term, which corresponds to 
YHLii n i ~ 1i) 2 + M/2 = 0, then applying (|105|l we find that 



/o=n C qi ) {-If'ZlM) ■ F q (0). 

i=l ^ qlJ 

The explicit form of F q (0) can be obtained from the following results 



(108) 



D„(0) = (-1)1^1^^ H ]Jp^(^), 

i=i 

1, Hcrmitc, 



G,(0) = 



HM= II (A*fc-A*i) 3 * m x 

1 <]<k<m 



2 M 2 £ 2=i«? 5 Hermite, 
1, Laguerre. 



(109) 
(110) 
(111) 



Now let us put together what we have learned about the asymptotic behavior oi lN,m,q(n) to describe the asymptotic 
behavior of HN,m,q(ti) when q € N m . To this end, we substitute (|109|l into (|108jl . then (|108jl into ()104|l . and finally 
substitute (|104|l and (jHEJ into (JHIJ, to obtain 



n 



N,m,q 



(M) = 7V Ef =I «(«-!) + JJ (/i 



1 G(2 9i , 

= 1 v ^ ' l<j<fc<m 

m 



l + O 



AT. 



(112) 
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Using the explicit forms for G q (0) and H q (Q), given by (|110f> and l|lll[) respectively, one can easily verify that for 
both the Hermite and Laguerre cases we have the following identity 



JJ e - 29 ^Rc { S^ (AiiK7 rf . Gq{0)Hq{0) . hQ 



i=l 



=n^( t , t )p- q H^)2 2q <-^^- 2q > n (pk-^ q ". (H3) 

l<j <k<m 



i=l 



Inserting the identity (|113|) into (|1 121) we finally obtain 

m m ^2 ( | -i \ 



(27T) 



1=1 

x n iMfc-^ii 

l<j<fc<m 



1 + 



(114) 



which does indeed recover (|21|l . 

We emphasize that the derivation we have presented for l|114|) is entirely rigorous for q £ N m for any m £ N, thus 
verifying the legitimacy of the log-gas procedure used ir£ for the Hermite and Laguerre cases, for such q. 
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